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O^- Abstract 

■^ . We explore theoretically electromagnetically-induced transparency (EIT) in a superconducting quantum 

circuit (SQC). The system is a persistent-current flux qubit biased in a A configuration. Previously [Phys. 
rj^ Rev. Lett. 93, 087003 (2004)], we showed that an ideally-prepared EIT system provides a sensitive means to 

Q ■ probe decoherence. Here, we extend this work by exploring the effects of imperfect dark-state preparation 

o : 

and specific decoherence mechanisms (population loss via tunneling, pure dephasing, and incoherent pop- 
►^ ^ ulation exchange). We find an initial, rapid population loss from the A system for an imperfectly prepared 

f— >) ' dark state. This is followed by a slower population loss due to both the detuning of the microwave fields 

t:j- , from the EIT resonance and the existing decoherence mechanisms. We find analytic expressions for the 

o 

\G ' slow loss rate, with coefficients that depend on the particular decoherence mechanisms, thereby providing a 

P ' 

means to probe, identify, and quantify various sources of decoherence with EIT. We go beyond the rotating 

2 ■ wave approximation to consider how strong microwave fields can induce additional off -resonant transitions 

'^ ! in the SQC, and we show how these effects can be mitigated by compensation of the resulting AC Stark 



shifts. 

PACS numbers: 82.25.-j, 03.67.-a, 42.50.Gy 



I. INTRODUCTION 



Superconducting quantum circuits (SQCs) based on Josephson junctions (JJs) exhibit macro- 
scopic quantum-coherent phenomena yj]. These circuits exhibit quantized flux or charge 
states states, depending on their fabrication parameters. The quantized states are analogous to 
the quantized internal (hyperfine and Zeeman) levels in an atom, and the SQCs thus behave 
like "artificial atoms." Spectroscopy IS Ij, |4|, |5|], Rabi oscillations and Ramsey interferome- 
try llajv, sLsl [lOl lil|,ll2|], cavity quantum electrodynamics llJ, Il4|], and Stiickelberg oscilla- 
tions lll5L 11611 are examples of quantum-mechanical behavior first realized in atomic systems that 
have also been recently demonstrated with SQCs. 

We recently leveraged the atom-SQC analogy to propose electromagnetically induced trans- 
parency (EIT) lIlTL llSh in superconducting circuits [ 119 0. EIT has attracted much attention in 



^ 



2311 and nonlinear 



atomic systems in the context of slow light [20], quantum memory 112 IL 
optics 1I24I1 . EIT occurs in so-called "A-systems" comprising two meta-stable states, each coupled 
via resonant electromagnetic fields to a third, excited state. For particular initial states called "dark 
states," the absorption on both transitions is suppressed due to destructive quantum interference, 
thus making the atom transparent to the applied fields. Though EIT is often studied in the context 
of the behavior of a weak 'probe' field in the presence of a stronger 'pump' field, we focus on the 
case where the two fields have comparable amplitude. In Ref. lll9ll . we analyzed a superconduct- 
ing persistent-current qubit biased such that it exhibited a A -configuration: two meta-stable states 
(the qubit) and a third, shorter-lived state (the readout state). We showed that EIT provides a non- 
destructive means to confirm preparation of an arbitrary superposition state of the qubit. Moreover, 
we showed that the proposed EIT scheme can sensitively probe the qubit decoherence rate using 
a method analogous to the proposal in Ref. |25] for atomic systems. This method compliments 
other available techniques of probing decoherence such as spin echo [26] and Rabi oscillation de- 
cay ['27']. Because the EIT method requires no manipulation of the qubit during the probing, it 
offers unique advantages in this regard. In addition to our EIT work, several groups have consid- 
ered the use of "dark states" in SQCs comprising a A-co nfigura tion to implement adiabatic passage 



and its application to quantum information processing 1128 
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In the present work, we extend and augment our analysis in Ref. lll9ll with realistic effects 
which arise in SQCs due to the presence of additional quantized levels (beyond the three-level 
"A-system" model). These effects have qualitatively unique signatures in an EIT experiment, and 



this work provides a tool for identifying their origin. This allows a more complete understand- 
ing of the full level-structure of the SQC system, and it further clarifies the necessary criteria for 
the experimental observation of EIT. The present work carries the spirit of previous investiga- 
tions in which additional degrees of freedom (beyond two-level models) were required to explain 



quantitatively experimental Rabi oscillations in SQCs. Examples of t 



tunneling across the barrier Bill , diagonal dipole matrix elements 



lese works include resonant 



- >2ll , and coupling to additional 



33ll . Just as EIT is sensitive to 



degrees of freedom outside the SQC, such as micro-resonators 11271 

decoherence, it will be similarly sensitive to effects beyond the idealized three-level model. 

The effects we investigate arise primarily from differences between SQCs and the atomic sys- 
tems considered in much of the literature. First, while damping of the excited level is provided 
naturally by spontaneous emission in atoms, in SQCs, this decay is 'manufactured' by resonant 
biasing across the tunnel barrier followed by fast measurement with a SQUID. This process must 
be considered in more detail to assure this decay is indeed analogous to spontaneous decay in 
atoms. Second, the transitions are at microwave rather than optical frequencies, whereas the Rabi 
frequency coupling rates and dephasing rates tend to be faster than in atomic systems. This means 
that various couplings in the system can be more comparable to the level spacings and thus the 
rotating wave approximation (RWA) is often not as valid as in atomic systems. Third, the level 
structure itself is quite different. In particular, there is typically some degree of dipole-like cou- 
pling between all pairs of levels in the system, because selection rules allow all possible transi- 
tions 1I34I1 . Fourth, in SQCs, there is the possibility of direct resonant tunneling across the barrier, 



a feature which is absent in atomic systems. 

This paper is organized in the followin g rn anner. In Section HIl we introduce the proposed 
system, a persistent-current (PC) qubit 11351 13611 . We discuss the conditions under which the PC 
qubit exhibits a A-configuration amongst its multi-level energy band structure that is conducive 
for an EIT demonstration. We then present the Hamiltonian and density matrix approaches to 
analyze the system dynamics. In Section [nil we use the Hamiltonian approach to give useful 
analytic approximations to the full system, which allow us to investigate EIT in a reduced three- 
level system. We explore effects of population and phase mismatch between the prepared initial 
state and the desired dark state (as defined by the applied fields), and the effect of detuning the 
applied fields from their resonances. We also consider the SQUID measurement rate and its effect 
on the effective decay and frequency of the excited 'read-out' level. In Section |Wl we use the 
density matrix approach to include pure dephasing and incoherent population loss and exchange. 



generalizing our previous results in Ref. lll9ll . We explore the effect on EIT in the presence of 
coherent and incoherent tunneling processes. Generally, one must make the EIT 'preparation 
rate' (proportional to the microwave field intensities) faster than the decoherence rate in order to 
observe EIT. In Section |V1 we go beyong the rotating- wave approximation (RWA) to examine the 
important issue of microwave field-induced off-resonant transitions in the spirit of previous work 
on two-level systems BTI] . We conclude that off-resonant transitions cause frequency shifts and 
losses which depend on the coupling field intensities. Unlike decoherence and detuning, these 
transitions generally manifest themselves as the field intensities are increased. The high-field 
frequency shifts are analogous to the AC-stark shifts observed in atomic systems. We show how 
the off-resonant effects can be mitigated and, in some cases, eliminated by compensating for the 
frequency shifts. 

II. THE PC-QUBIT 

A. The PC-qubit model 

The persistent-current (PC) qubit is a superconductive loop interrupted by two Josephson 



junctions of equa 



(Fig. Eta)) 



ml 



size and a third junction scaled smaller in area by the factor 0.5 < a < 1 
. Its dynamics are described by the Hamiltonian 
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+ Ej [2 + a — 2 cos (y9p cos (y9m — ttcos(27r/ + 2iprn)] , (1) 



in which C is the capacitance of the larger junctions, (^p „ = (y^i ±y92)/2, ^i is the gauge-invariant 
phase across the larger junctions i = {1,2}, Ej = Ic^q/2tx is the Josephson coupling energy, Ic 
is the critical current of the larger junctions, and / is the magnetic flux through the loop in units 
of the flux quantum $o Bal . 

The qubit potential energy (the second term in Tipc) forms forms a 2D periodic double well 
potential, a one-dimensional slice through which is shown in Fig.lTJb). Each well corresponds 
to a distinct classical state of the electric current, i.e., left or right circulation through the loop, 
with a net magnetization of that is discemable using a dc SQUID [36]. The relative depth of 
the two wells can be adjusted by detuning the flux bias to either side of the symmetric point 
/ = 1/2. The potential wells exhibit quantized energy levels corresponding to the quantum states 
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FIG. 1: A pc-qubit with a dc SQUID measuring device (a) A pc-qubit, a superconducting loop with 
two Josephson junctions of equal dimension and the third scaled by a factor a, as shown in the inner 
loop. The outer loop is a dc-SQUID that is used to measure the magnetic moment of the qubit. (b) A 
representation of the potential energy of the pc-qubit as a function of /, the magnetic flux in the loop in 
units of the flux quantum ^q. The qubit potential can range from an asymmetric double well biased to 
the right to a symmetric double well and to the an symmetric double well biased to the left for / ranging 
from < 1/2, = 1/2, and > 1/2 respectively, (c) One-dimensional double-well potential and energy-level 
diagram for / = 0.502, in which case we have a three-level system in the left-hand well. States |1) and 
1 2) are meta-stable, while |3) will have significant loss via resonant tunneling to |4) ((T34). The right- 
hand well states undergo fast damping (r4,r5,r6) via the SQUID measurement and intrawell relaxation 
to lower states. Coupling between the our three-levels is induced by two resonant microwave fields with 
Rabi frequencies ili3 and 1^23. forming the A-system. The qubit parameters we use in calculations are 
UJ2 — ^1 = (27r) 27.8 GHz and L03 — 1^2 = (2vr) 27 GHz, with matrix elements Xij for (i, j) = (1,2), (2,3), 
and (1,3) set to -0.0145, -0.0371, -0.0263, respectively. 



of the macroscopic circulating current 



uy 



40, 



4111 . with the number of levels on each side 



determined by the depth and frequencies of the wells. In this basis the Hamiltonian can be written 



(2) 
(3) 
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— Tto + /ttunneli 
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'T-tunnel 


= n''^aij\i){j\ 




i,j¥=i 



(4) 



We note here that two points of view may be taken when discussing the system described by 
the Hamiltonian in Eq. @. In one picture, the diabatic states (diagonal matrix elements) of the 
qubit are the uncoupled single- well states of classical circulating current, and these diabatic states 
are coupled through the tunneling terms (off-diagonal matrix elements). In a second picture, the 
Hamiltonian is diagonalized, resulting in eigenenergies and eigenstates of the double- well poten- 
tial. Although the perspectives differ, these two pictures will, of course, lead to identical results; 
only the interpretation differs. Throughout the paper, we primarily describe the dynamics in terms 
of diabatic (single- well) states coupled through the tunneling barrier of the double- well potential 
and driven by harmonic excitation. Exceptions, where they exist, will be clearly noted. 

The three-level A structure to implement EIT is then provided by the left-hand meta-stable 
states |1), |2) and the fast decaying level |3) shown in Fig. [He). Each of these levels are taken 
to have a finite loss rate r| ', due to resonant tunneling to a right-well state (at aij) followed by 
relaxation of the right-hand well state Tj (which is a sum of population relaxation to lower levels 
and damping induced by a fast SQUID measurement of the circulation current of right-hand well 
states). In particular, we desire a fast decay rate r3 , which is achieved by resonantly biasing |3) 
and 1 4) and a fast SQUID measurement (~ 1-10 ns), as analyzed in Section ITlIDl Conversely, we 
desire states |1) and |2) to be long-lived and the tunneling (T25 will cause loss and decoherence, 
which is analyzed in Section HVDK aif; is negligible by comparison). Rough estimates of the 
interwell loss rates for resonant-tunneling are F-^ ^ (1 ms)~^, F2 (1 yus)^^, and F3 ^ (1 ns)^^ 
and the off -resonant biasing of states 1 1) — |6) and |2) — |5) will significantly decrease these rates. 
In addition, states |2) and |3) can have intrawell relaxation rates T^^i, T-i^2, r2^i (not shown in 
the diagram). Under similar bias conditions the rate F3_,i + F3^2 ~ (25 /is)^^ (experimentally 
measured in ll40ll ) is much slower than F3 . Also note that, F2^i, another source of decoherence 
of the meta-stable states, will be less than T^^i + F3^2- 



These quantized levels may be coupled using microwave radiation. An applied radiation field fx 
can be described in terms of an amplitude, frequency and phase: A/^ = g^ cos^tu^t + cp^). We find 
the resulting matrix elements for level transitions (the Rabi frequencies) by treating A/^ as a small 
perturbation in the cos(27r/ + 2^9^) term in Eq. dH). We write it as sin(27r/ + 2ipm) sin(27rA/^), 
which can be approximated as sin(27r/ + 2(/)m)(27rA/^), leading to a Rabi frequecy uf^' = 
{2n)gfj_aEjXij/h, where Xij = {i\ sin(27r/ + 2Lpm)\j)- The elements Xij we calculate for our 
proposed parameters are listed in the caption of Fig. [He). In EIT, we address the SQC with two 
microwave fields A fa = Qa cos{uJat + (f)a), with cUa ~ ^3 — ^1, and Af^ = Qb cos{uObt + 0b), with 
ujb ^ LOz — ^2- The microwave induced Hamiltonian can then be written as: 

H.-wave = I E E (fi!f e-^(-^*+^^) + c.c) \^) {j\ (5) 

where i,j runs over the states and n runs over the two fields a, b. We emphasize that the above 
approximation is a perturbative approach valid only for small driving amplitudes. In the strongly 
driven limit, the approximation breaks down, preventing the Rabi frequency from growing without 
bound libj]. 

Microwave excitation is used to establish the population of meta-stable states (such as |1) and 
|2)) \ia photon-assisted tunneling. In this scheme, the population of a meta-stable state is driven 
via a resonant radiation field into a read-out state {i.e. |3)) which quickly tunnels to a measurement 
state (i.e. |4)). This state has opposing current circulation with a unique flux signature that can be 



measured using a DC-SQUID 113611 . In the present scheme, two fields are simultaneously applied 
(resonant with |1) ^^ |3) and |2) <-h> |3); see Fig.[Hc)), and EIT is manifested by a suppression of 
the photon-assisted tunneling due to quantum interference between the two excitation processes. 

B. Evolution model 

It is convenient to calculate dynamics from the above Hamiltonian terms in an interaction pic- 
ture which transforms away the diagonal energies Hq Q. In this frame the total Hamiltonian is 
then the sum of Q and ^: 






Note that the exponential arguments involve sums of microwave frequencies u^ and level splittings 
Ui — LUj. When these nearly cancel the state is said to be near-resonant and the coupling is strong. 
However, for most of the terms, this cancellation does not occur, and the term rotates its phase 
rapidly on the scale of frequencies of interest (^L , cTij and Tj). Such terms are neglected in the 
Rotating Wave Approximation (RWA). 

We can also include incoherent losses from the levels, F^, by introducing an additional non- 
Hermitian par t of the Hamiltonian TYreiax = — ^^ J2i(^i/'^)\^){A- This is often done in quantum 
optics 1I42L I43I1 to include non-Hermitian decay of radiatively decaying levels. We then describe 
the system by a wavefunction |\^) = J2i Ci{t)\i), with the initial population normalized to unity 
X]j l^i(0)P = 1 (this can decay in time due to the non-Hermitian loss). The evolution of the |^) 
is the governed from Schrodinger's equation: 

^^^|^) = (7i: + 7i:,eiax)|^). (7) 

Besides giving the coherent dynamics, this Schrodinger equation correctly predicts the population 
relaxation of level \i) at Fj and also gives the correct dephasings of coherences between \i) and 
other states at half this rate Fj/2. 

When necessary, we use a density matrix approach to include incoherent processes. For ex- 
ample, pure dephasing of a coherence between the two meta-stable states |1) and |2) goes beyond 
the Hamiltonian approach ©. Similarly, incoherent feeding of levels (such as population into 
|1) from interwell relaxation F2^i) goes beyond this description. The density matrix is written 
P ~ J2ii PijV) iJ I ' where the population in the levels are given by the diagonal terms pu and corre- 
spond to |cjp in the wavefunction description, while the off-diagonal terms pij correspond to CiC* 
and describe coherences between levels. The evolution of the density matrix is given by 

d ^ ~ ~ 

^^ olP = [^ + ^relax, P] + 'Cp. (8) 

The first term reproduces the part already predicted by the Schrodinger equation (jT]), while the 
super-operator C, the Lindbladian [l42l], accounts for other incoherent processes. For pure dephas- 
ing of the \i) <-> |j) coherence, 7jj, we introduce a term Cij^ij = —712- For a population relaxation 
from \j) -^ \i), Tj^i, we introduce Cjj^u = +Fj_^j. The associated population loss from \j) and 
decoherences are already included through Eq.|7](via a term ~ih{T j^i/2)\i) {i\). 

Throughout the paper, we consider the model in a number of distinct cases. In each, we include 
three levels |1), |2), |3), coupled by two microwave fields A/^, A/^ making up our A system (see 

8 



Fig. [He)). We then selectively include additional levels, such as the |4), |5), and |e) to isolate 
the contributions of each of them. Numerical results were obtained with a fourth-order Runge- 
Kutta algorithm ll44ll solving Eq. ^. In it we do not make any RWA assumptions a priori, but 
instead introduce some cut-off frequency ujrwa- We examine the phase factors of each term in the 
evolution and set to zero ones with phases rotating faster than ujrwa- 

We compare our numerical results with approximate analytic solutions in many cases. When 
possible, we use the Schrodinger equation dTJ to obtain simpler analytic results, though the full 
density matrix approach is used when dephasing and interwell relaxation are considered (Sec- 
tions |lVAl-|lVCl)- In the analytic results we normally make an additional transformation |\^) -^ 
\'^), H+Hic\a.x -^ 7Y, defined by transformations of each level frequency {q} -^ {q} = {cjC*''**}, 
where the 5i are chosen in to eliminate time-dependent exponential phase factors in ® (they are 
usually detunings, that is, frequency mismatches between the microwave frequencies and the cor- 
responding transitions). Detuning from two-photon resonance, decoherence, and additional levels 
are all seen to destroy the perfect transparency of EIT and cause slow exponential loss of the 
population. We will obtain expressions for the loss rate Rl in these cases. 

III. ELECTROMAGNETICALLY INDUCED TRANSPARENCY IN A SQC 

A. Ideal EIT in a A configuration 

We first consider the 'ideal' case in which the three levels in the left well (see Fig. [He)) are 
well isolated from direct tunneling to other levels, states |1) and |2) are perfectly stable, and |3) 
quickly decays at some fast rate Fg . This decay is in reality due to resonant coupling of |4) (0-34) 
and subsequent SQUID measurement r4, but we will see in Section ITlI DI how one can derive F^ ' 



3 



in terms of these underlying processes. 

We apply fields with nearly resonant frequencies cua = to^ ~ cui + A13 and cuh = uj^ — uj2 + 
A23 (see Fig. [He)), where the A13, A23 are small detunings. All other couplings are sufficiently 
detuned to safely eliminate them under the RWA. In this case the transformations to eliminate 
phase rotating terms are given by 61 = 0, ^2 = A13 — A23, ^3 = A13. The Hamiltonian, written 
in matrix notation in a basis {|1), |2), |3)} is 



^-2 









0* 





-2A2 


^^23 


fil3 


^23 


-zr« - 2Ai3 



(9) 



where A2 = A13 — A23 is the detuning from two-photon resonance. Here we have dropped the 
a, b labels, ^13 = ^ij and ^^23 = ^23 (^^^ Eq. Q) as there is no ambiguity. The open system 
loss of 1 3) due to tunneling Fg is assumed to dominate incoherent population exchange due to 
intra-well relaxation, allowing a Schrodinger evolution analysis Q. 

First consider the resonant case A13 = A23 = 0. A qubit initially in the ground state |1) can 
be prepared in a superposition state |\l'init) = ci|l) + C2|2) by temporarily driving it with a field 
resonant with the |1) ^^ |2) transition. Applying only one field ^13 (^23) field then allows the 
population of a state |1) (|2)) to be read out through a transition to state |3) followed by a rapid 
escape to the right well. In this case, the superposition is destroyed by the absorption of a photon. 

However, from O and ® it follows that if we simultaneously apply both fields and the SQC 
is in the dark state 

|vl^.) = ^|l)-^|2), (10) 



(where f2 = A/|fii3P + 1^23^) then dV]) predicts ci = C2 = C3 = 0. For this particular state, the 
two absorption processes, |1) ^^ |2) and |1) ^^ 1 3), have equal and opposite probability amplitude 
and thus cancel by quantum interference. As a result, no excitation into |3), and thus no tunneling 
to the right well, will be observed. Note that \^d) constrains both the relative intensity and phase 
of the light fields. Any other (non-zero) values for the relative amplitudes in the two states 1 1), |2) 
will lead to a coupling into |3) and subsequent loss. 

An alternative interpretation is obtained by examining the eigensystem of the Hamiltonian Q. 
The dark state I^E'd) has eigenvalue zero. The other two eigenstates are linear combinations of the 
excited state |3) and the combination of the stable states orthogonal to |^d): |^a) = (^i3|l) + 
^2312))/^, called the absorbing state (a "bright state"). The system (I^^a), |3)} acts effectively 
as a two-level system coupled by Vl. The eigenvalues corresponding to the two eignestates are 
(— ir3 ± \l AVP — Fg )/4 and the imaginary parts of these eigenvalues give the loss rates of 
these states. In the limit Vl ^ F3 , these rates are both F3 /4 and one observes damped Rabi 
oscillations. In the limit Vl <^ F3 , there is an eigenstate ~ \^ a) with a slower damping rate 
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FIG. 2: (Color online) Suppression of tunneling due to EIT for various ideal wavefunctions (a) The 

populations of the states as a function of time in the presence of applied fields Qi^ = ^223 = (2vr) 150 MHz 
and tunneling rate Fg = {2tt) 130 MHz = 1/1.2 ns for the initial state pu = P22 = 0.5 , pi2 = —0.5 
(the dark state). The dotted (red) curve shows pn, the dashed (blue) p22 and the thin solid (green) p^^. The 
total population (sum of the three) is the thick solid (black) curve. There is a slow exponential decay of the 
population due to the dephasing rate 712 = (27r)l MHz. (b) The population evolutions (same convention) 
for the initial state pu = P22 = 0.5 , pi2 = 0.5 (the absorbing state), (c) The population evolutions for an 
initial state |1) (which is an equal superposition of the dark and absorbing states). 

Figure Eta) shows an example of the lack of tunneling in the presence of applied fields ^13 = 
r223 for the corresponding dark state |^init) = (|1) — |2))/"\/2 (i.e. pu = P22 = 0.5 , pi2 = —0.5). 
One sees only a barely perceptible population P33 and a very slow loss of the pn and p22- This 
is due to a pure dephasing of the state coherence, which we take to be 712 = (27r) 1 MHz. The 
effect of this dephasing is a small exponential loss at a rate we label R^ , which is discussed and 
derived in \\% and reviewed in Section lTV Bl Otherwise the populations remain pn = P22 ~ 0.5. 
EIT thus provides a means to confirm, without disturbing the system, that one had indeed prepared 
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the qubit in a particular desired state of the SQC, preserving its quantum coherence. 

By contrast, Fig. |2tb) shows the large loss induced when one applies these same fields to the 
absorbing state, i.e., the state with the same populations but n out of phase: |\E'init) = (|1) + 
|2))/v^. In Fig.|3b) we see that there is a large population in the |3) and the entire population 
has tunneled to the right well within about 10 ns. Note that here we are in the intermediate 
regime fi ~ Fg so we get oscillations with period ~ i7 strongly damped at ~ 1^3/2. This is 
the completely analogous to the tunneling which occurs with a single applied field in a two-level 
scheme. 

A general state can be decomposed into dark and absorbing state components. Fig. Etc) shows 
a case where the initial population is purely in |1) and the same fields are applied. Here the initial 
state can be written |\E'init) = |1) = (I^d) + \^a))/V^- Half of the population (the component 
in the absorbing state) is coupled out over the 10 ns time scale while the dark state component 
remains. In terms of level populations pn, P22, approximately 1/4 of the population is coherently 
coupled from |1) to |2). 

B. EIT with imperfect state preparation 

One of the useful aspects of EIT is the extremely sensitive manner in which it can measure the 
amplitude and phase of superpositions in the SQC. When the prepared state has a slightly differ- 
ent phase or population ratio than the state we intend to prepare, EIT could be used to measure 
these deviations. Such imperfect preparation could arise, for example, due to imperfections in the 
preparation pulse. 

FigEta), shows the populations loss when preparing a state with various initial state population 
ratios and applying fields rtn = ^23- We again introduce a small dephasing 712 = (27r)l MHz. 
The inset shows the population at 200 ns, well after the initial transient losses have occured. This 
data can be understood using the dark/absorbing basis discussed above. The modulus square of 
the overlap of the initial state and the dark state (^Dl^init) gives the population remaining after 
the fast initial loss of the absorbing component. Postulating that the slower loss (due to dephasing 
or other effects) is exponential with some rate Rl, the population after the fast initial loss is: 

|(^init|^D)|V^^*. (11) 

Thus, detecting the fast initial decay of the population can indicate the mismatch between the fields 
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FIG. 3: (Color online) Imperfect state preparation(a) The time evolution of the population in the left 
well as a function of initial state of the form |^init) = \/Pi\^) — \/(l — Pi)|2). The uppermost curve is for 
the dark state pi = 0.5. Successively lower curves are for pi =0.6, 0.7, 0.8, 0.9 and 1.0. There is a sharp 
initial decay when there are deviations from the dark state. Inset: The population in the left well at 200 ns 
versus pi. (b) The population decay out of the left well as a function of the initial phase of the prepared 
state, l^init) = (|1) - e*^|2))/\/2. 9 = (top curve) and the other curves are vr/S, 2tt/5, 3tt/5, in/b, and 
TT. We see full decay for = tt, the absorbing state \^a) = (|1) + |2))/\/2. Inset: The population in the 
left well at 200 ns as a function of 9. 

and the prepared state population. 

Phase mismatch, or unwanted 2; -rotation, in the qubit preparation shows similar behavior. Fig- 
ure|3tb) shows the population decay from the left well for the state (|1) — e*^|2))/-\/2. The upper 
most line is decay due to the perfect state, while the lower lines indicate the decay for varying 
value of 9. This example indicates that the dark state, is more sensitive to phase mismatch than 
population mismatch. 

C. EIT detuned from resonance 

Because EIT is a coherent effect, it only occurs in a narrow range of frequencies near the two- 
photon resonance. The width of the EIT feature is generally determined by the field intensities, 
and can be made narrower than the broad resonances of the individual one-photon transitions 
(|1) ^^ |3) and |2) ^-> |3)), which are determined by the fast decay rate of |3) Fg . 

Figure|4ta) shows the results of simulations with the same parameters as Fig. 12 a), but with the 
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FIG. 4: (Color online) EIT in the presence of detuning, (a) Numerical calculation of the total population 
in time when A13 = MHz, and at various A23 (top to bottom curve) (27r) 0, 10, 20, 30, 40 MHz. 
For A23 = 0, the decay is due to pure dephasing, while the decay is sharper when A23 7^ 0. Inset: The 
population in the left well at 200 ns as a function of the two-photon detuning A2 = A13— A23, with A13 = 
(circles). The solid curve shows the prediction (fT3b . The diamonds show the case A13 = (27r)l GHz, 
varying A2 about the two-photon resonance, (b) The population decay as a function one-photon detuning 
at two-photon resonance, A13 = A23 = A for (top to bottom) {2tt) 3, 2, 1, 0.5, GHz. As the one-photon 
detuning increases, the effective coupling of the fields to the transition reduces, hence reducing the rate of 
decay. Inset: The population in the left well at 200 ns as a function of the detuning for the dark state (circles). 
For comparison we also show the population remaining for the absorbing state |^init) = (|1) + |2))/\/2 at 
the same detunings (diamonds). 



detuning A23 varied. The curves show exponential loss occuring at various rates. We can analyze 
the results with the Hamiltonian ® and the corresponding Schrodinger equation Q. We first 
adiabatically eliminate ll45ll the excited level by setting C3 = and obtain 



C3 = (fii3Ci + VI23C2) 



2A 



13 



iT 



(t) 



4A?3 



r 



(t)2 



(12) 



This expression is valid for times long compared to the initial transient time MinKFg ^^ , A]]^^}. 
Note that for the dark state (ITOt the amplitude C3 vanishes. Plugging this expression back into 
the equations for ci , C2 then gives a 2 x 2 matrix evolution equation, which can be easily solved 
by finding for its eigenvalues and eigenvectors. For A2 = the eigenvectors are simply the dark 
l^u) and absorbing {"^a) states of Section JIIIAL with eigenvalues Xd = and A^ = —Q'^{r^^' — 
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2zAi3)/2(4Af3 + r3 ), respectively. The absorbing component population is damped out at 
— 2i?e{Ayi}. In many cases, we are interested in the regime close to the one-photon resonance 
Ai3 < r^*^ for which this reduces to Q'^/T^^^ as in Section ing 

With a non-zero two-photon detuning A2, this process is complicated an additional phase evo- 
lution term. Plugging Q into Schrodinger's equation dTJ) (in the frame defined before ^) gives a 
term C2 = • • ■ + « A2C2, which acts to drive the phase of the SQC out of the dark state and competes 
with damping of the absorbing component. Solving for the eigensystem in this case we see that, in 
the limit of small two-photon detuning (IA2I ^ |i?e{A^|}), the eigenvalue corresponding to the 
dark component has a non-zero negative real component, leading to a population loss rate: 

in |2in |2 a2f(*) 
-2Re{Xo} ^ Rf^' = ^\M^^ (13) 

The prediction P = exp[— (i?}^^ + R\^ ^ )t] is plotted in the inset of Fig.|4ta) and is seen to agree 
well with the numerical results (where i?)^^^ was determined by the numerically calculated loss 
for A23 = 0). Eq. (fT3t shows how the field strength, via Q,'^ in the denominator, determines the 
frequency width of the EIT feature. 

The above analysis indicates that it is only the two-photon detuning A2 which effects the 
relative phase of |1) and |2) and therefore effects the dark state. EIT will occur in the pres- 
ence of a large one-photon detuning A13 and the inset of Figure Ha) shows such a case with 
Ai3 = (2n) 1 GHz and A23 varied about the two-photon resonance. The presence of a trans- 
parency peak is still clear. The important difference is, because of the large one-photon detuning 
Ai3 ^ Fg , the damping of the absorbing state — 2i?e{A^} is substantially reduced and so both 
the dephasing R^^^^' and detuning R]^ ^' loss are reduced. The analytic model (IT51) is not valid 
for large one-photon detunings A13 where the strong damping assumption —2Re{\A} ^ IA2I 
does not hold. As a result, one sees non-exponential decay in the large one-photon detuning cases 
(upper curves of Fig.|4tb). 

Figure lltb) shows simulations at two-photon resonance A13 = A23, varying the one-photon 
detuning A13. The population decay is much slower as the detuning gets larger. For comparison, 
we also plot the decay for initial state equal to the absorbing state |\l'init) = (|1) + |2))/-\/2. We 
note that the analytic model for loss of the dark state (fT3t is invalid for large one-photon detunings, 
where the absorbing state is not completely damped. 
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D. Effective A-system via tunneling and measurement 

Thus far we have considered the system to be a three level system with the excited level |3) 
subject to a fast population decay Fg . Underlying this decay are actually two processes: the fast 
resonant tunneling to a near degenerate level in the right hand well ((T34) followed by interwell 
relaxation and possibly a strong measurement of the population in |4) (r4); see Fig. [He). We show 
here how the picture of a three-level system with a strong damping of |3) (the Hamiltonian Q) is 
most valid when (T34 < F4 but actually has a larger range of validity than one might expect. We 
derive an expression for F3 and also see how the tunneling slightly shifts ^^3. 

To do this we consider the Schrodinger evolution of the full four-level system Hamiltonian 
(with the same frame transformation as Eq. ® and 64 = A13 + ^34, where ^34 = 104 — lu^): 
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(14) 



We have used the notation A\.J to distinguish it from A13 which includes the frequency shift of cus 
induced by |4). 

To recover our three-level picture, we note that when F4 ^ (T34 we can adiabatically elim- 
inate level |4) to obtain C4 = —2030-34/ [2(^34 + A{^) — iF4]. Plugging this result back into 
the equation for C3 reveals that our system can be reduced to a three-level system as in ^ with 
F« = 4|a34pF4/[F2 + 4(^34 + Ag))^] and A13 = AJ? + Mas4\%,/[rl + 4(^34 + AJ^V]. Al- 
tematively, when F4 <^ 0-34 we would expect the tunneling to induce a splitting of |3) and |4) into 
two superposition eigenstates (split by 20-34). 

We carried out several numerical simulations of the four- level density matrix equations dH) for 
this system, considering first the resonant case (^34 = 0). In them we used CT34 = (27r)150 MHz = 
(1.2 nsy^, fields Qrs = (27r)120 MHz, ^23 = (27r)150 MHz, the corresponding dark initial 
state pii = 0.61, P22 = 0.39, pi2 = —^ p\\P22 i}-^- full coherence), and a dephasing rate 712 = 
(27r)2 MHz. Figure |5ta) shows the populations P33,p44 versus time for cases Fj^ = 1 ns = 
((27r)159 MHz)-i (thin, blue curves) and 25 ns = ((27r)6 MHz)-^ (thick, red curves). Fig.gtb) 
shows the total population remaining versus time. In the fast measurement case, 1 ns, P33 and 
P44 are seen to track each other, and we see the exponential decay of the population as in the 
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FIG. 5: (Color online) Consequences of the measurement state characteristics, (a) The thinner curves 
(blue) show populations p^^ (solid curves) and p^s (dashed) as a function of time for a fast read-out T J^ = 
1 ns < a^^ (with scale on the right side). After initial transient period, the two values reach quasi-steady 
state values, which undergo slow exponential decay. Conversely the thick (red online) curves show a slow 
measurement case FJ^ = 25 ns » a^^ (scale on the left). In this case the populations do not reach a 
quasi-steady over the time scale plotted, (b) The total population remaining versus time for the same two 
cases, (c) The population remaining at 100 ns for varying measurement rates (dots) and compared to the 
prediction of a three-level system with F3 . As the measurement gets slower, this prediction slightly and 
increasingly underestimates the actual population which should be observed, (d) Numerical results (dots) 
and three-level model predictions (solid curves) now letting ^34 vary, for the cases F J^ = 1 ns (black, lower 
curve) and 10 ns (upper, red curve). 

previous examples. For the slower measurement case, 25 ns, we see p^^ and P44 still track each 
other, but here there is a large excitation P33 (note the different scale), no fast transients in P33 to a 
quasi- steady state value, and non-exponential population decay, all indications of the breakdown 
of BIT. 

To check the validity of the effective three-level model, in Fig.|5tc) we compare its predictions 
for the populations remaining after 100 ns (solid curve) with predictions of the full four- level 
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model (dots). The agreement is excellent for the 1/Ta, < 5 ns, and is still in rough agreement even 
up to 50 ns. The adiabatic elimination procedure appears to be valid well beyond the expected 
regime r4 ^ ^ 0-34^. The breakdown of EIT in the 25 ns case is due to T^^' becoming too large (see 
Section ITVB I ), rather than a breakdown of the effective three-level model. 

In Fig. |3td) we show the population remaining for when ^34 is varied (for both l/V/^ = 1 ns 
(black curves) and 10 ns (red curves)). We have kept the microwave fields on bare resonance 
A ^3^ = A23 = and accounted for the predicted frequency shifts of u^ in the three-level model. 
As (534 becomes comparable to r4 the tunneling rate Fg ' goes down as predicted. 

In summary, we find that so long as F4 > (734 the simple three-level provides an excellent model 
and even when F4 ~ (134 or somewhat larger, this model unexpectedly gives very good predictions 
of the behavior. However, one must be careful of the strong dependence of F3 ' on F4, as this may 
severely effect the necessary conditions for EIT (which are discussed in Section ITVB I) . Through 
F4 the SQUID measurement rate can thus have a large influence on the EIT. When the detuning 
(534 becomes comparable to F4 the tunneling rate is reduced as expected and one must account for 
the frequency shift of uj^. For the remainder of the paper, we will not explicitly include |4) in the 
calculations, but assume some Fg ~ 1 ns and that the frequency shift is already included in the 



definitions of A 13, A 



23- 



IV. EIT IN THE PRESENCE OF DECOHERENCE, INCOHERENT POPULATION EX- 
CHANGE, AND QUBIT TUNNELING 

An outstanding, important issue in the eventual application of SQCs to quantum computing is 
the identification and suppression of sources of decoherence and unwanted dynamics of the qubit 
states (|1) and |2)). In particular, decoherence from pure dephasing ll46ll . intrawell relaxation, 
interwell resonant tunneling [31], and coupling to microscopic degrees of freedom in the junction 
ll27h have all been proposed as potential hurdles in successfully isolating a coherent two-level 
system for use as a qubit. Use of phase sensitive methods such as EIT could be a fruitful path 
for exploring and differentiating the contributions of these various decoherence processes to qubit 
dynamics. To learn how EIT is effected by decoherence, we use the density matrix approach here 
to include pure dephasing, population loss, incoherent population exchange, and resonant coupling 
to the right well. We find a minimum microwave coupling strength necessary for the observation 
of EIT in the presence of decoherence and see it contributes small exponential loss (as seen in the 
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numerical results above). We derive analytic expressions for the loss rates, which are proportional 
to the decoherence processes but with coefficients which depend on the nature of the process. 



These results are a generalization of the results for pure dephasing previously published lll9ll 



EIT is a unique tool to probe decoherence which compliments the previously explored tech- 
niques of spin echo ['26'] and Rabi oscillation decay f2r\, and this section will demonstrate some 
of its advantages. First, in the limit of small decoherence, the system is left completely undis- 
turbed by the probe. Second, the dependence of the loss rates on relative field strengths can be 
used to determine the nature of the decoherence process. Third, the populations of the qubit states 
|1) |2) do not need to me manipulated in the process (with n pulses, etc.). Besides the advantage 
of simplicity, this latter point also leaves open the possibility of decoherence rates which have a 
non-trivial dependence on the relative state populations. 

A. Density matrix approach 

To carry out this analysis, we must go beyond the Schrodinger approach, and introduce the 



corresponding density matrix evolution ©, also referred to as the Bloch equations ll42ll . We work 
in the three-level case and transform to the frame defined above Eq. Q and obtain: 



Pll = -ri Pll + r2^lP22 - 2^13P31 + 2^13Pl3, 

P22 = — (r2 + r2^l)p22 — ^^23P32 + ^^23P23; 

P33 = — r3 P33 + -r2]^3P3i — -f2i3Pi3 + -fi23P32 — ^^23P23! 

P12 = -(712 H ^ JP12- 2^^13P32+ 2^^23Pl3, 

Pl3 = ^-^-^Pl3 + 2^t3(Pll - P33) + 2^23Pl2, 

P23 = ^ P23+ 2^23(P22-P33) + 2^13P21- (15) 

For simplicity, we have supposed the resonant case A13 = A23 = and ignored inter-well re- 
laxation of 1 3), which is dominated by Fg . The remaining three elements equations are de- 
termined by p*j = pji. The most important new piece here is the decoherence rate of P12: 
7i2 + (ri + Fg + F2^i)/2. This will decohere the dark state and lead to the small losses in 
the population. 
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To proceed, we adiabatically eliminate the excited state coherencess pia, P23 as they are strongly 
mped by Fg . In these equations we can ignore F-^ , r2 ^ F; 
then plug the results back into the remaining equations to obtain: 



damped by Fg . In these equations we can ignore F-[ , F2 ^ F3 as well as P33 ^ pn, p22- We 



(i) l^'isl /''l3*'23 
Pll = -ri Vll + F2^lP22 f^Pll - ( -^p P12 + C.C), 

, (J) , ''23 / * '13* '23 \ 

P22 = -(r^ + r2^l)P22 :^ P22 - I -1^ P12 + C.c), 

P12 = -(712 + ^ )pi2 ^Pi2 oT^^ii + ^22) (16) 



2 >-- 2F«"^^ 2Fr 



(«) 



We note here a strong damping of the coherence provided by the fields r2^/2F3 ' . This damping acts 
to drive the system into the dark state. In the limit that the decoherence terms (712, T\\ Fg , F2--*i) 
vanish, there is a steady state solution consisting of perfect coherence: pn = |i723p/f2^, P22 = 
1^13 P/i7^, P12 = —^13^23/^^- However, the decoherence terms drive the system out of the dark 
state, causing excitation P33. One sees the ratio of the decoherence rate compared with the EIT 
preparation rate Vl^ /V^^' determines the degree to which the coherence deviates from the perfect 
dark state value. 



B. Measuring dephasing with EIT 

We first show this comes into play for the pure dephasing, which is expected to be the case 
in many practical implementations and was analyzed previously in [19]. Figure |6t a) (solid, red 
curve) shows the excited state population P33 when 712 = (27r)l MHz and we apply the fields 
Vti^ = ^23 = (27r)150 MHz to the dark state |^init) = |*d) = (|1) - |2))/v^. One sees a 
small (note the scale in Fig.|6ta)) but finite excitation. In particular, there is an initial fast transient 
behavior to some plateau value (over a time scale determined by Minjfi^/Fg , F3 }), followed by 
a slow exponential decay. The general behavior of an initial transient rise, Fig.|6ta), was seen over 
a wide parameter regime. This quasi-steady state excitation of |3) is the origin of the exponential 
losses at rate w^^'^' in the previous simulations. Fig.l^b) shows the exponential decay for several 
different dephasing rates 712. The inset of Fig. |6tb) plots the populations remaining at 200 ns, 
which is seen to approach unity as 712 -^ 0. 

The loss rate can be quantified by considering the 3x3 evolution matrix for pn, P22, and pi2 
defined by (IT6b . Looking at the eigenvalue corresponding to the smallest loss rate, and expanding 
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FIG. 6: (Color online) EIT loss due to pure dephasing (a) The population P33 versus time in the presence 
of a pure dephasing 712 = (27r)l MHz (solid, red curve, scale on left) and 712 = (27r)20 MHz (dashed, 
blue curve, scale on right). In the slow dephasing case P33 quickly reaches a plateau value ^3^^^^ then 
undergoes a slow exponential decay. For the faster dephasing, the exponential decay time is similar to 
the time required to reach the maximum value. In each case, the initial state is pu = P22 = 0.5 and 
p^2 = -0.5, ili3 = i723 = (27r)150 MHz and T^*^ = (27r)130 MHz (b) The population decay with varying 
dephasing rates (top curve to bottom curve) 712 = {2tt) 1, 2, 3, 4, and 5 MHz. Inset: The population in 



the left well at 200 ns versus 712 (dots). The solid curve shows the analytic prediction exp{—R}J^^ t) based 
on Eq.[n] demonstrating how the population loss can be used as a probe of 712. (c) Population remaining 
at 100 ns versus detuning A13 (keeping A23 = 0) for two different field intensities. In these simulations 
we used the initial conditions pn = 0.61, P22 = 0.39, pu = —^ P\\P22, the fields $^13 = 0.8(723, 
Fg = (27r)159 MHz, 712 = (27r)l MHz. The solid curves are the analytic prediction described in the 
text for $723 = (27r)100 MHz (blue, upper curve), and 30 MHz (black, lower curve). The dots show the 
numerical results. 

to first order in 712 gives: 



n 



{712) 



10 |2|0 |2 
I "13 1 I "23 1 



^L - 2712 ' '"'y (17) 

By measuring this decay constant experimentally, one can use Wh to extract 712. Note that the 
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two rates are simply related by a constant or order unity, determined by the relative strength of the 
two microwave field couplings (in the example in Fig. |6l the constant is 2^2^3^223/^'^ = 0.5). A 
glance at Fig.l^b) reveals how choosing the observation time t ~ 7^2^ will give the best sensitivity 
in the measurement. 

The inset of Fig. I^b) shows the prediction ( flTl in comparison with the analytic results and 
we see good agreement. The adiabatic elimination and the expansion for small 712 require 
27i2r^*V^^ < 1 and 712 < r^*^ This ratio is 0.07 for 712 = {2n) 5 MHz. For higher de- 
phasing rates, the dephasing rate competes with the preparation rate, and the first order expansion 
in 7i2 becomes less valid. Such a case is seen in Fig.|6ta) (dashed, blue curve) where we plot P33 
for a case with 712 = (27r) 20 MHz. The exponential decay occurs with a time scale comparable 
to the transient time to reach the quasi-steady state plateau. In such cases, 712 can only be esti- 
mated from the tunneling rate by more detailed modeling of the underlying Bloch equations (ITSb . 
We note that the microwave field intensity can be adjusted to control fi^ and to bring us into a 
regime where (flTt is valid. The breakdown of this inequality occured in the slow measurement 
time (25 ns) case plotted in Fig.|5ta,b), for which r3 ^ 4|cr34p/r4 became quite large (~ GHz). 

When detuning and dephasing are both present but sufficiently small, the two effects add lin- 
early. In Fig.|6tc) we show the population remaining versus detuning A13, with 712 = (27r)l MHz 
(with A23 = 0) both in a strong field (27i2r(*)/fi2 = 0.OI6) and weak field (27121^/^^ = 0.18) 
case. The prediction P = exp{—{R£^^' + R]^ ^')t) holds for the stronger field but overestimates 
the loss for the weaker field. 

We conclude from the above calculations that there are two basic conditions that have be satis- 
fied for a reliable measurement of the decoherence in the system. First, the decoherence rate 712 
should be much smaller than the loss rate Fg , which holds in systems of interest. Second, we 
must be able to apply sufficiently strong fields that the preparation rate fi^/Fg dominates 712. 

C. EIT with incoherent population loss and exchange 

When the decoherence of pi2 occurs due to population loss and exchange instead of dephasing, 
the effect on the EIT is much the same, with the 712 simply replaced by the total decoherence rate. 
However, because these processes involve additional changes in the populations, the population 
loss rate, which we use to diagnose the decoherence rate, will be different. 

Referring back to evolution matrix (O, we again find the eigenvalues to determine the loss 
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FIG. 7: (Color online) EIT loss with dephasing, open loss, and closed loss The population loss after 
100 ns for several types of decoherence. The parameters are as in Fig.|5Jc) but with il23 = (27r)150 MHz. 
The curves show the analytic predictions ( fTTl and (IT^ and the dots show numerical results. The lowest 
curve is for purely open loss Fg , the middle curve for pure dephasing 712, and the top curve shows pure 
population exchange r2^i. The horizontal axis shows the corresponding decoherence rate for each case: 
r^*V2, 712, F2-.1/2, respectively. 

rate of the dark state. Expanding to first order in F^ \ r2 , r2-^i, respectively, we find: 
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It is interesting to note that the coefficient will depend in different ways on the relative intensities 
of the two fields depending on the origin of the decoherence. Fig.|71shows the population loss after 
100 ns for different kinds of loss, each plotted versus the total decoherence rate of pi2. The open 
system loss r2 is greater than the pure dephasing case because there is direct population loss on 
top of the absorption into P33 due to decay out of the dark state. The closed system loss (intrawell 
relaxation r2^i) is seen to be smaller than pure dephasing, however, (ITSt shows that this could be 
greater or smaller depending on the relative values of ^13 and ^23- 

While we assumed in this discussion that the intra- well relaxation of |3) was negligible (r3^i + 
r3^2 ^ Tg '), the present discussion is easily generalized when this can not be assumed. In this 
case the population loss rate is the same as before, multiplied by a factor reflecting the proportion 
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(t) 



of atoms in |3) which actually tunnels to the right well: r3 /(Fa^i + r3^2)- This reflects the fact 
that excitation of |3) due to decoherence will only be registered as population loss upon tunneling 
to the right well. 



D. Resonant tunneling loss out of left well 



Just as the decay of |3) is the result of tunneling followed by decay of |4), interwell tunneling 
of 1 2), which is conceivably a leading order effect in the decoherence, involves tunneling to a near 
resonant level |5) (see Fig.[TJc)). Here we consider this effect in detail to find the conditions where 
an effective damping rate r2 can be used, and also explore conditions where the dynamics are 
more complicated. 

To explore this issue, we have performed numerical simulations for the system 
{|1), |2), |3), |5)} where level |2) is detuned from |5) by ^25 = cjs — uj2. The Hamiltonian (with 
the transformation ^5 = A2 + S25) is: 
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20-25 



2(725 -^15 - 2(A2 + fe)). 

Mathematically similar energy level structures have been considered in the context of atomic sys- 
tems [47]. 

In the following we take 0-25 = (27r)5 MHz and set the pure dephasing 712 = to isolate the 



contribution from the presently considered effect. Analogous to Section IIII Dl when Fs ^ (725 
one can easily reduce the system to an effective-three level system with an additional loss rate 
F2 = 4|a25pF5/(F5 + 4^25). The small (black) dots in Fig. Eta) present the population as a 
function of the detuning A13 (keeping A23 = 0) after 100 ns of evolution for a case with T^^ = 
2 ns ~ (15 0^25)^ (and ^25 = 0). The results are in good agreement with the prediction one obtains 
from the loss rate Eq. (fTSt with this predicted tunneling rate Fg (solid curve). Fig.[8tb) shows 
the computed loss (dots) for resonant fields (A13 = 0) compared with the loss expected from the 
calculated Fg (solid curve) as a function of Fg^^ for two different field strengths. One sees the 
estimate is good for T^^ < 5 ns. The inset shows P55 versus time when T^^ = 5 ns. One sees it 
quickly reaches a quasi-steady state plateau, then undergoes an exponential decay. 
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FIG. 8: (Color online) Loss due to resonant tunneling to right well. In the simulations we assume a 
(resonant) tunneling rate 0-25 = (27r)5 MHz, with fields i7i3 = 0.8i723 and the corresponding dark state 
pii = 0.61, P22 = 0.39, p\2 = —^ P\\P22- (a) The population remaining at 100 ns of fields applied 
with strength Q.2Z = 50 MHz, versus the detuning A13 (keeping A23 = 0), but varing the relaxation time 
r^^. For the small solid dots (black) F^"^ = 2 ns. The solid curve shows the analytic prediction based 
on the effective loss rate V^ described in the text. The open (blue) dots show a case F^^ = 15 ns, in 
which case a splitting appears in the resonance, contrary to the analytic prediction (dashed curve). The 
large (red) dots show the case F^^ = 80 ns, for which the splitting becomes more pronounced and the 
loss rate quite small, while the r2 model predicts complete loss of the population, (b) The population 
remaining at 100 ns at the two photon resonance (A13 = A23 = 0) versus F5" . The solid (black) dots 
show the case f^is = (27r)150 MHz and the open (red) dots show ^13 = (27r)50 MHz. They roughly agree 
with each other and the Fg model (solid curve) for F^"^ ^ al^^ = 32 ns. However, for larger F^"^ the 
loss becomes slower. The inset shows the population P55 for the a fast (5 ns, solid curve) and slow (80 ns, 
dashed) relaxation times F^"^ (note the different scales). The fast case looks analogous to decoherence (see 
Fig.|6ta)), while oscillations occur in the slow case, (c) The population remaining versus the level detuning 
^25 in the case V^ =8 ns (solid, black dots) and 80 ns (open, red dots). The solid and dashed curves show 
the F2 model predictions. 25 



For larger F^^, Fig.[8tb) shows the loss begins to decrease in contrast to the analytic estimate. 
The red curve in Fig.[8ta) shows the population remaining versus detuning A 13 in a case in for 
F5 = 80 ns. Besides the Fg model incorrectly predicting complete loss of the population after 
100 ns, in the numerical results there is the clear appearance of double-peaked structure. This can 
be understood from the coupling (725 giving rise to two eigenstates (|2) ± |5))/v^ split by 2(725, 
each of which gives rise to a distinct BIT resonance. The initial state (no population in |5)) is a 
superposition of these eigenstates and we get oscillations of the population between |2) and |5). 
The inset in Fig.[8tb) (dashed curve) shows these oscillations in P55 for F^^ = 80 ns. Because of 
the weak damping, the oscillations persist and the quasi-steady state is not reached during the time 
scale plotted. The dotted curve and open dots in Fig.[8tb) show an intermediate case T^^ = 15 ns 
where the double peak structure is just becoming apparent, and the analytic estimate has begun to 
break down. 

In Fig. lite) we address the case where the tunneling levels |2) and |5) can be slightly off- 
resonant {§25 7^ 0). The filled dots show the population remaining for F^^ = 8 ns versus ^25. The 
solid curve shows the Fg model estimate, which correctly accounts for the slower tunneling rate 
as we move off resonance. The red dots show the same for F^^^ = 80 ns. In this limit, the analytic 
estimate expression severely overestimates the loss for S25 < F5 but as we move off resonance, 
the coherent tunneling plays less of a role and the effective tunneling decay rate model becomes 
more accurate. 

In summary, tunneling of our lower states will be a source of loss in EIT. The behavior will 
depend qualitatively on the relative strength of the coherent coupling and the loss rate of the 
additional quantum level and so can provide us with information about these quantities. In the 
limit where the loss rate dominates, we see how the it reduces to an open system loss of |2) where 
as in the other limit we see a qualitative signature (the splitting of the resonance) of coherent 
coupling to another level. 

V. EIT WITH RADIATION CROSS-TALK 

To now, we have considered how EIT is affected by decoherence and tunneling to other levels. 
Another important consideration to include is that all levels are dipole coupled and so in principle 
coupled by the microwave fields. The RWA allows us to neglect the majority of the couplings as 
the dynamics are dominated by couplings which are near resonant. For example, one need not 
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consider the coupling of field b on the |1) ^^ |3) transition or field a on |2) ^-> |3). However, in 
SQCs, the relative scale of the Rabi frequencies (~ 100 MHz) to the level spacings (~ GHz) is 
somewhat larger than in typical atomic systems. Thus, it is important to know the magnitude and 
type of effects that these "cross" couplings can have. Here we consider, separately, a case with 
cross coupling within the three level system, and a case where fields couple to an additional excited 
level. In general, we find these effects can be characterized analytically in terms of additional loss 
rates and AC Stark shifts. If one neglected these effects, there are configurations where one may 
mistakenly attribute a loss rate to a dephasing when in fact it is due to off-resonant field coupling. 
While we found it was beneficial to turn the microwave coupling strengths Vl^ up to overcome 
decoherence and detuning, we will see how this can increase the importance of these cross-talk 
effects. We will also see how proper understanding of the effects can allow us be mitigate them 
by properly compensating for the Stark shifts. To isolate these cross-talk effects, we will set other 
losses and dephasings to zero in the following. 

A. Radiation cross-talk in a three-level system 

In the configuration proposed here and in lll9ll . the qubit states 1 1) and |2) are the first two levels 
of a slightly anharmonic potential (the left well), while the excited level |3) is the third such level. 
Therefore, the level spacing u^ — uj2 is only slightly different than uj2 — uj\ (with the parameters 
proposed the spacings are ~30 GHz and the difference is 0.7 GHz). As a result, the field ^^23 is 
only 0.7 GHz detuned from |1) ^^ |2) (see Fig.|9ta)). A rough estimate of this effect was noted 
in 1I19I1 . Here we present an analytic treatment which shows it causes an AC Stark shift which 
depends on the relative dipole coupling strengths, field intensities, and level spacings. Thus the 
EIT resonance position can be a function of field amplitudes used, which can be compensated by 
adjusting the field frequencies. 

We consider the Hamiltonian ® but do not invoke the RWA with respect to terms rotating by 
the mismatch frequency between the |2) ^^ |3) and |1) ^^ |2) transitions 5 = {uj2—^\) — [,^-i—^'2)- 






(3*n*,e''' ni. 



pn23e-''' -2A2 n*2s 

fii3 ^23 -2Ai3 - ^r^*^ 



(20) 



where (3 = Xi2/x2'i is the ratio of dipole moments between the additional off-resonant transition 

27 




5 10 15 20 25 30 35 

time(ns) 





.-R-^J^ 




S0.9 


..^^r-*' / 


V *^^ 


•i-H 


J» S' • / 


\ * 




.•'5' / 


\* '« 


3 0.8 


J' •/ 


\ '* 


a 


- / / 


\ • 


o 


• / 


\ 


^0,7 


1 


\ • 








C^ 




\ 


-4— > 


1 


\ 


firiO.e 


1 ml 


\ • 



■10 -5 



-5 5 

^J2n (MHz) 



10 



FIG. 9: (Color online) Cross talk in a ladder system, (a) Schematic of the dominant cross-talk term: field 
h (resonant with |2) <-^ |3)) also couples |1) <-> |2), detuned by 0.7 GHz. (b) This induces fast oscillations 
of the ground state populations pn and P22 (and a slow overall drift) as shown here for the initial state 
pii = 0.61, P22 = 0.39 with full coherence and O13 = 0.8 O23 = (27r)120MHzandrf^ = (27r)159 MHz. 
(c) Population remaining at 50 ns for i723 = (27r)50 MHz (small, blue dots), 100 MHz (open, green dots), 
and 150 MHz (large, red dots), with the curves showing the analytic predictions based on the AC Stark 
shifts described in the text. 

and the intended resonant transition for the field h. In the case we are considering /5 = —2.55 
though it should be emphasized that these ratios are strong functions of the parameters and can 
vary by an order of magnitude. 

We performed a numerical propagation of the density matrix equations ^ for this Hamiltonian 
for the resonant case A13 = A23 = and plot the evolution of pn and P22 in Fig.|9tb). We see the 
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extra coupling gives rise to a small amplitude, rapid oscillations of both quantities. Considering 
a toy two-level model with only the off-resonant coupling present predicts population oscillations 
of period {2ii)/5 and amplitude ~ 1/9^231/25, in agreement with the numerical results. The small 
deviations of pn, P22 from their dark state values gives rise to absorption into |3) and thus loss. In 
the toy model, an off -resonant coupling can be accounted for as an AC Stark shift. In particular, 
uji and UJ2 are predicted to shift by ± ~ \l3VL2-i\'^ / A5 , respectively. This results in an effective shift 
of the two-photon detuning A2 which can be compensated. 

Stated in terms of our exponential loss language, the loss rate R)^ (O is still valid but the 

(12) 
two-photon detuning A2 should be replaced by A2 + A^^% where 



^(12) ^ ^^^^2^ ^ 

AC 25 ^ ^ 

In Fig. |9tc) we plot the population remaining versus A13 (keeping A23 = 0) for three different 
values of field intensities. The solid curves show the theoretical prediction based on the predicted 
AC Stark shift. They are in good agreement (the overestimate of loss at the lowest intensities 
with some detuning is due to the damping of the absorbing that being too weak to efficiently keep 
the SQC in the dark state). Importantly, if one adjusts the field frequencies, one can completely 
avoid loss due to the cross-talk coupling. Strictly speaking, there is a small loss if |2) decays at 
some small rate Fg , however, this loss is much smaller than the loss already predicted from the 
associated decoherence (ITSt . 

B. Effect of off-resonant radiation coupling to an additional excited level 

The last important situation we consider is that of coupling to an additional excited level |e), 
coupled off-resonantly to states |1) and |2) via the two applied fields a and h, thus forming a 
"double-A" system 11481] . as diagrammed in Fig. flOl' a). As we will see, the extra coupling gives 
rise to AC Stark shifts in much the same way as we saw in Section fV^ In addition, because |e) 
(unlike |2)) is quickly decaying, the coupling gives rise to some population loss even when the AC 
Stark shift is compensated. 

Dropping the RWA with respect to terms coupling to level |e) and defining 5-ie = oje — oJ-i^ the 
Hamiltonian (using a frame 4 = A13 + 5^e) is 
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FIG. 10: (Color online) Coupling to an additional excited level, (a) Schematic of off -resonant microwave 
coupling of each of states |1) and |2) to an additional level |e) above the barrier, (b) The population 
remaining at 100 ns for the same initial state and relative field strengths as in Fig. ^ We show the cases 
^23 = (27r)60 MHz (small sohd, blue), 100 MHz (open, green), and 150 MHz (large solid, red). The solid 
curves show the loss and AC Stark shift predicted in the text ( l^5t . To isolate and clearly show the effect we 
have set /3 = (from d^ ) and used 63e = (27r)1.5 GHz, instead of the (27r)10 GHz we predict for our 
proposed parameters. We use xu = 0.0054 and a;24 = -0.0437 and Tg = T^*^ = (27r)159 MHz. 
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with Pi = Xii/xis. We have assumed some large open loss channel Tg. 

In the case where only one of the couplings is present (/3i = or /92 = 0), the effect is simple 
to calculate. When /3i = 0, one can consider the Schrodinger evolution Cg from Eq. (|22ll and 
adiabatically eliminate Ce to obtain 



■ /^2^23C2 



(23) 



3e 



where we have assumed S^e ^ A13. Substituting this back into the equation for C2 yields: 



C2 = --fi23C3 - 'iA2C2 - 1/^2^ 



231 



2i53e + Te 



C2 



(24) 



revealing that the extra coupling gives rise to a Stark shift and population decay of |2). In the large 
detuning limit (d^e ^ Te), the Stark shift is |/32^23p/453e and we have an effective loss rate from 
|2), Fg = 1/92^^23 pre/4(5|g. Analogous results occur when [32 = 0, leading to a Stark shift of 
the two-photon resonance of opposite sign. The population loss rates will in turn contribute to the 
decoherence and cause exponential loss from the EIT as discussed in Section HVl 

When both are present, the shifts and loss rates are not simply the sum of the two separate 
contribution, due to interferences between them. A case where summing the two contributions 
clearly does not work is Pi = P2 (equal dipole ratios). In this case (and only this case) the dark 
state C2/C1 = —^13/^23 is also completely decoupled from |4). Thus no population loss or Stark 
shift is induced in this case. To obtain an expression in the general case, we follow the following 
procedure. We adiabatically eliminate C3 and Cg from the Schrodinger equation from Eq. (E^ . 
obtain the 2x2 evolution matrix for ci and C2, then solve for the eignvalues of this matrix. One 
of the values has a large imaginary part (which reduced to f^^/Fg in the limit Pi = P2 = 0) 
and corresponds to the absorbing state. The other has a small imaginary part which vanishes 
when Pi = P2 = and corresponds to the dark state. We investigated the imaginary part of 
this eigenvalue in the limit that fi^/Fg :^ A13, A23 and PiVe/S^e ^ 1 and then minimized this 
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expression with respect to the two-photon detuning A2 (when A23 = 0) to obtain the Stark shift. 
The loss rate and shift of the resonance obtained were: 

in |2in |2 n2F 

p(e) _ |"13| |"23| " J-e , n n \2 

Rl - — ^, — Isi-^f'-^'\ ' 

Aji = -^^|^(A - /52)(A|fil3p + /52|fi23r). (25) 

These expressions reduce to the simpler cases above (/3i = or /92 = 0) and also disappear when 
Pi = P2- Interestingly, the relative sign of the dipole moments plays an important role. For 
example if Pi = —P2 the loss rate is actually twice what one would expect from the sum of the 
individual couplings. In this case, the dark state \E'd for |3) is actually the absorbing state for |4). 

In Fig. \W[ h) we plot the population remaining versus detuning A 13 for three different field 
intensities in the system we have been considering, but now considering the coupling to the addi- 
tional level |e). The dots show numerical propagation of the density matrix equations correspond- 
ing to Eq. (l22b . Note that to isolate the effect studied at present, we have ignored the cross-talk 
considered in Section fV Al (bv setting X12 = 0) and set the pure dephasing 712 = 0. The solid 
curves then show the analytic estimate based on the Stark shifts and loss rates (l25t . One sees 
excellent agreement. 

It should be noted that this analysis should be able to account for the effect of multiple excited 
levels Cj by simply summing their contributions ^jR^ and X]j ^ac- Because of the large 
frequency differences between each successive level, coherent interference between contributions 
from different \ej) will not occur. 

VI. CONCLUSION 

We have described in detail a proposal for demonstrating a quantum optical effect, FIT, in a 
SQC. In this context, FIT will manifest itself as the suppression of photon-induced tunneling from 
stable states |1), |2) through some read-out state |3), due to quantum mechanical interference for 
two paths of excitation. This provides a method of unambiguously demonstrating phase coher- 
ence in these systems. We have provided a thorough and mostly analytic treatment of FIT in 
the presence of complicating effects due to decoherence and multiple levels in SQCs, which will 
be important in guiding experimental implementation and observation of FIT and other quantum 
intereference effects. 
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We analyzed in detail first the basic considerations of EIT such as imperfect dark state prepa- 
ration, and one- and two-photon detuning and determined the expected experimental signatures. 
Under appropriate conditions, we obtained an expression for the total population as a function 
of time Eq. dTTT) . which describes a fast loss of the absorbing component, followed by a small 
exponential loss of the system. Over shorter times, EIT thus provides a method to confirm the 
successful preparation (and coherence) of the particular dark state defined by the microwave fields 
applied. For longer times, the observed loss rate will be function of both the detuning from two- 
photon resonance as well as decoherence effects. We also discussed the important issue of how 
the measurement state |4) plays a role in the decay of the read-out state |3) and saw how the bi- 
asing condition of these levels and the SQUID measurement rate can have a large effect on the 
parameters of the effective three-level system. 

We then discussed in detail how decoherence due to dephasing of the qubit coherence, inco- 
herent population loss or exchange, and tunneling of levels through the barrier effects the loss 
rate. Measuring these loss rates can then be a powerful tool which sensitively probes these various 
processes. We obtained the coefficients for the loss rates, which depend differently on the field 
strengths, depending on the underlying decoherence processes. For the case of primarily coherent 
resonant tunneling, we found that that the EIT will exhibit a qualitatively different double-peaked 
structure. Probing these effects with EIT can aid in understanding and minimizing decoherence 
and give information about the full multi-level structure of the SQC. A potentially interesting fu- 
ture investigation is to learn the signature from coupling to other quantum degrees of freedom, 
such as the microresonators postulated in ETIl . 

Finally we have found that the microwave fields themselves can cause additional loss rates 
and AC Stark shifts of the EIT resonance which must be accounted for when one uses stronger 
field strengths. Importantly, we found that, these effects can becomes more pronounced with ^2^, 
meaning there will be some intermediate field strength which balances these considerations with 
the decoherence and detuning effects. Also, we showed how these effects could be mitigated by 
proper compensation of the Stark shifts. 
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